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ERROR IN INDUCED ELECTRIC FIELD MODES UNDER STIX COIL 
RESULTING FROM ASSUMPTION OF ORTHOGONALITY OF 
NATURAL PLASMA MODES 
by Richard R. Woollett 
Lewis Research Center 

SUMMARY 

When the Bessel function representing the induced electric field in a plasma under a 
Stix coil is expanded in terms of the natural plasma modes, the coefficients of the expan- 
sion may be approximated by assuming that the eigenfunctions are orthogonal. However, 
the eigenfunctions are not orthogonal. When the orthogonality assumption is not used, 
there results an infinite system of equations involving an infinite set of unknowns. Kantor- 
ovich developed a technique that can be adapted to this problem. The technique produces 
both a deficient and an excessive value of the coefficient. Because the correct value lies 
between these two limits, the greatest error resulting from the orthogonality assumption 
can be obtained. In the cases tested it was found that the errors involved in determining 
the coefficients using the orthogonality assumption were less than 0. 1 percent for the 
fundamental mode. 


INTRODUCTION 

When the efficiency of adding energy to a plasma wave by means of a Stix coil (ref. 1) 
was calculated, a Bessel function in the analysis was expanded into an infinite series of 
eigenfunctions (refs. 1 and 2) as follows: 


- E a m J l<V R > 

m=l 


(Symbols are defined in appendix A. ) The coefficients of the infinite series were then de- 



termined by the usual Fourier series technique, which requires that the eigenfunctions be 
orthogonal. This straightforward approach yielded an easily evaluated, closed-form ex- 
pression for the coefficients: 


a 


m 



RJ 1 (yR)J 1 (y m R)dR 



Rjl(u m R)dR 


The eigenfunctions, however, are not orthogonal; consequently, the cross product terms 
are not identically zero, 


J RJ i( J ' m R ) J i( l ' n R ) dR ^ 0 

and the preceding expression for the a m 's is incorrect. 

There exist techniques for solving problems involving an infinite set of nonorthogonal 
eigenfunctions, such as the method of Green's functions (ref. 3) and Schmidt's orthogonal- 
ization process (ref. 4). These require expansion of the original nonorthogonal functions 
into polynomials which form an orthogonal set. However, these techniques result in ex- 
pressions which are difficult to evaluate and compare with the coefficients in the original 
series, that is, with the a m 's in 


oo 


L 

m=l 


a m J l(% R ) 


However, by proceeding in a manner similar to that used with the Fourier-Bessel expan- 
sion, a direct comparison can be made. In this latter technique, the cross product terms 
do not vanish, and, instead of the simple expression for a m presented previously, there 
results an infinite set of simultaneous equations, 


f 0 1 RJi( l .R)J 1 (^R)dR=2^a m ^ 'rj^RJJ 


(^ n R)dR 


m=l 


for which it is not possible to obtain an exact solution. However, by an iteration method 
described in reference 5, it is possible to calculate upper and lower limits of the coeffi- 
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cients. If the iteration process is extended, these two limits can be brought as close to- 
gether as desired. In this manner, a direct comparison can be made between the or- 
thogonal and nonorthogonal coefficients. This report calculates these upper and lower 
values of the coefficients for several plasma conditions near ion cyclotron resonance. In 
addition, these limits on the coefficients are compared with the values of the coefficients 
obtained from the orthogonal expansion used in references 1 and 2. 


THEORY 

The physical problem in reference 1 concerned the coupling of radiofrequency power 
into a plasma by generating ion cyclotron waves. A cylindrical plasma of infinite length, 
confined by a steady longitudinal magnetic field and surrounded by a vacuum, is capable 
of exhibiting many natural modes of oscillation. Among these are oscillations identified 
as hydromagnetic or ion cyclotron waves (ref. 6). The radial variation of the electric 
field in the plasma for the mth mode of oscillation is proportional to a cylindrical Bessel 
function of the first order and first kind, J^(f m R). The eigenvalue v , which is the 
radial wave number of the mth natural mode of oscillation, is determined from the plas- 
ma dispersion relation and the boundary condition at the vacuum-plasma interface (ref . 6). 
The axial wave number k of the mth natural mode is simultaneously found. 

If a disturbance in the plasma column is driven by a current sheet (Stix coil) periodic 
in both time and space, the radial variation of the electric field in the plasma is propor- 
tional to the Jj(i'R) where v is the radial wave number associated with the driving elec- 
tric field and is determined from the dispersion relation and the axial wavelength of the 
current sheet. 

In reference 1 the forced electric field consequently was written as 


OO 

EccJ 1 ( v r)= £ Vit^R) (1) 

m=l 

the electric field of the mth natural mode being proportional to 

a m J l^m R ) 

It is necessary to consider this expansion because the efficiency of power transfer was 
expressed for each natural mode separately. If a m can be evaluated for nonorthogonal 
J l(^m R ) and com P ared with the coefficient obtained when the orthogonal assumption is 
made, the maximum error of the coefficients can be determined. Because the electric 
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fields of the natural modes are proportional to the a m 's, the relative errors of the a m ’s 
are also the relative errors of the field modes. To calculate the a m 's, multiply equa- 
tion (1) by RJ^ n R) and integrate with respect to R over the range 0 < R < Rp where 
Rj is the radius of the plasma. Then, 



1 

RJ 1 (i/ n R)J 1 (yR)dR = 



R Jl(^ n R )Jl(%R)dH + a n 



RJ^(i/ n R)dR 


( 2 ) 


If the Jj(^ n R)’s are assumed orthogonal, all integrals of the 


r R l 

T RJ^RJJ j(^ m R)dR n £ 


m 


form vanish. However, these integrals can be expressed as 


/ R 1 


RJ 1 (^ n R)J 1 (i/ m R)dR = 


R l J l(v m R i) J i(^ n R l) 


2 2 

V - v 
m n 






Since, in general, the eigenfunctions are nonorthogonal, the right side cannot be expected 
to vanish for all n and m combinations. If the eigenfunctions were orthogonal, either 
Jj(^ n Rj) would be zero or the term ^ n Rjjj 0 (^ n R^)/j^(^ n R^)j would be constant for all n. 
All cross products would then vanish. However, J^z^R^) does not vanish, and 

is actually a slowly varying function of i^Rj which approaches 
a finite limit. This latter condition is responsible for the small error calculated herein. 
The nth equation of the set represented by equation (2) becomes, after integration, 


^n ^ ^ ” ^nm^m^nm + a n"^nn 
m=l 


where 


n = JiC^R^J^^Ri) 

c n 


2 2 
" ^n 


v n R l 


J o< y n R l> 

J l^n R l> 


vRj 


J o^ R l) 

J^Rj) 


(3) 


(4) 
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J l(^ n R l) J l(V R l) 


nm 


2 2 

v - v 
m n 


v n R l 


Jjl' R,) 

o v n r f t* cr m V 


- v R* 
m 1 




m^n 


(5) 


_r| J^Rj) 

nn 0 9 

2 w 


v n R l 


J oW 


v12 


J l(^ n R l) 


- 2 v n R l 


J oW 

J l(^n R l) 


+ ("n R l>‘ 


( 6 ) 


and 


nm 


4 ln 

1.0 n 


= m 


t m 


Solving for a in equation (3) yields 


m=l nn nn 


In Kantorovich's notation (ref. 5), define 


C nm - (6 - 1) 

nm ' nm 


nm 


nn 


and 


b n" 


nn 


Hence, 


a = C a + b 
n Z-> nm m n 

m=l 


(7) 


If orthogonality is assumed, the sum term vanishes and the solution for a n is just b n , 
the Fourier-Bessel coefficient. The error is obtained by comparing the values of a n de- 
termined from equation (7) with the b n> 

In order to evaluate a , it is first desirable to determine the signs of the C ^ 

n nm 
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terms. Because all the subscripted v's that appear in A nm and A nn (eqs. (5) and (6)) 
are wave numbers of the natural modes, they satisfy the boundary condition (ref. 6) for 
undamped waves 


"n R l 


J oW 


= -k R, 

n 1 


K o< k n R l> 

Wi) 


( 8 ) 


where k = (/c^ - $^)V2 ^2 s, ^2 an( j \ and K i (k R 1 ) are modified Bessel func- 

n n ' ’ n ’ o n r 1 n 1 

tions of the second kind of order zero and one, respectively. Substituting the right side 
of equation (8) in equation (6) results in the following expression for A nn : 


_ r2 

nn 


2 K R i) 2 


r 


"N 

o 

< 

K (k R 1 ) 
k R. ° n 1 

Z KjkR.) 9 

+ 2k R. ° n 1 + (i'R) 2 > 


. K l< k n R l>. 

" 1 K^R,) “ J 


Because the modified Bessel functions are always positive, A nn is always positive and 
the sign of C nm will be opposite that of A^ m . To determine the signs of the A nm 's, 


consider equation (5). Because the term 


J o<"„ R l> 


is negative, the derivative 

J oWl 


^n R l 


J oW 


J l^n R l) 


^n R l 


k n R l 


Wi) 


K l( k n R i) 


v n R l 


1 n 1 i = 2 - ° . V n y . -(1^4) - {v n R t ) ° n 1 


d ( y n R 1 ) J l^n R l^ J l(^n R l) 

is also negative. Therefore, y n R i[j 0 ( t ' n R i)/Ji( l/ n R l)] is a rnonotonically decreasing 
function. Consequently, when m > n, that is, v > v , 


, W „ J 0 (y R i) 

'"■Wl* ^ J l<"m R l> 


> 0 


(9a) 


and 


2 2 . n 
"m" ^n > 0 


(9b) 
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Because the reverse inequality holds for both expressions when m < n, the ratio of the 
left side of these inequalities is always positive. The sign of A nm (eq. (5)) then depends 
only on the sign of Jj(i^ m R)Jj(^ n R). Because both 1 / and v must satisfy the boundary 
condition expressed in equation (8) and because both K o (k n R^) and K^k^R^) are greater 
than zero, J o (^ n Rp and J^(^ n Rj) must have opposite signs. It is shown in figure 1 that 
the nonparity of Jj and J Q occurs only in the bracketed regions between the zeros of 
the Bessel functions. 

Because the absolute magnitude of J Q /Jj varies from zero to infinity within these 
regions, there will be a v which falls within each of the bracketed regions. The value 
denoted by v ^ will be in the first region, in the second, etc. However, 

f> 0 n odd 

J lK R i) { 

1 n L< 0 n even 


This allows A nm to be expressed as 


A = (-1) 
nm ' 


n+m 


R l J l< 1 'm R l> J l< 

"n R l> 

, 2 - , 2 
m n 



J o<V R l> 


" hW ” J l< V R 1» 


( 10 ) 


Equation (10) explicitly exhibits the sign of A nm . Thus, C nm can be written as 

n+m+li 


c = (-i) n+m+1 |c 

nm ' I nm 


(id 


Particular care must be exercised when the signs of the coefficients vary. Rearrange 
equation (7) into the following form: 


N , 

N „ 

oo 

7 C nm a m + T 

nm m / j 

C a + b + C a 

nm m n Lj nm m 

m=l 

m=l 

m=N+l 


( 12 ) 


where N is a finite number, Y' indicates a sum of terms in which C nry _ > 0 (i. e. , 

nm 

n + m is odd), 2^ indicates a sum of terms in which C nm < 0 (i. e. , n + m is even), and 
the last term is a remainder term which will be approximated on page 8. 

An expression for a deficient solution can be synthesized from equation (12) by replac- 
ing the terms on the right side with expressions that are smaller. Similarly, an exces- 
sive solution may be obtained by substituting larger terms. The remainder term 
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£ 

m=N+l 


C a 
nm m 


is modified by selecting a number $ which is larger than every a m for every m>N + 1. 
The numerical technique for this selection is given in appendix B. The remainder term 
can be written as 


CO 


Yj C nm a m <^1? 
m=N+l 


where 

OO 

R ( n ) s y | C I (13) 

N I nml v ’ 

m=N+l 


If a n represents a deficient solution and a n an excessive one, write 


and 


N. 

N„ 

C mn 5 m + b „ -<**$ 




II 

M 

c_J + V 

nm m 

(n = 1, 2, . . 

N) 

(14) 

m=l 

m=l 





N 

N 

C nm a m + b n + '^n* 




E’ 

C a + Y" 
nm m / -> 

(n = 1, 2, . . 

., N) 

(15) 

m=l 

m=l 






Equations (14) and (15) represent 2N coupled equations in 2N unknowns. This set of equa- 
tions may be simplified somewhat by considering a n - a n and a R + a n as new variables. 
Hence, when equation (14) is subtracted from equation (15), there results 

N N 

K - - £' Cnm^m ' *m> + E" ' a m> + ^ 

m=l m=l 

However, because C nm in the Y," sum is always negative, 
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N 


S n - »„> = £ l C nml®m ' %> + <" = »- 2 • N > <«> 


m=l 


Similarly, when equations (14) and (15) are added, 


N .. 


<*„ + 5 n> = £ C nm< 5 m + V + £ + V + 2b n 

m=l m=l 


N 

£ 

m=l 


C (a + a ) + 2b 
nm' mm 7 n 


(n = 1, 2, . . . , N) 


(17) 


Equations (16) and (17) each consist of N equations and N unknowns, which may be 
solved by any one of many techniques. After values for a n - a n and a fi + a n are ob- 
tained, a n and a n are readily found. All that is required, consequently, is to increase 
N until and a n are nearly equal within desired limits. 

The solution of equations (16) and (17) depends on calculating values for the C nm ’s, 
which in turn depends on calculating many natural modes. However, the higher modes 
are very difficult to determine by normal techniques. These techniques generally break 
down because the solutions are near singular points. In addition, evaluating the C nm ’s 
requires calculating differences of nearly equal quantities (cf. , eq. (5)). This results in 
a considerable loss in computing accuracy. These difficulties may be circumvented, 
however, by expanding the Bessel functions in a manner that utilizes certain mathematical 
characteristics of the natural modes. To accomplish this, notice that the wave numbers 
for the natural modes approach the zeros of J Q (^ n R) as n increases without limit (ap- 
pendix C). It then becomes convenient to expand both J Q (^ n Rj) and J^(^ n Rj) about the 
nth zero of the zero-order Bessel function of the first kind. If j Qn represents the nth 
zero and t the expansion variable, it can be shown (see appendix D) that 


J oW * K + 7T- Fl(W 

^ J on/ 


and 




J on j 


on 


(18) 


(19) 
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where terms higher than the second power in t n have been neglected, 
which these are derived are convergent. The expression i^Rj J 0 (i/ R 
appears in C nm (see eqs. (5) and (6)) consequently may be written as 


The series from 
1 >/.I 1 (v rl R 1 )] that 


■'A 


- 4 i on t n(i 0 n + f) 


J oW _ 

J l (l 't R l ) 4J*„ - 4j on t n + 4t= - 


J on n 


( 20 ) 


because 


"n R l = U 


+ 


n 


(21) 


The variable t which appears in equations (18), (19), and (20) can be determined by 
utilizing the dispersion relation (ref. 6) 


r.4^2 »2 

£2 ot - 12 

2 OL+ (k 2 - fi 2 )_ 

(k 2 - S2 2 ) + (k 2 - O 2 ) 2 

o 2 

a + ( K 2 - O 2 ) 

- (k 2 - O 2 ) 


and the boundary condition 


i/Rj 


J 0 ( yR l) 

J 1 (t'R 1 ) 





( 22 ) 


(23) 


If equations (18) and (19) are substituted into equation (23) and use is made of the fact that 
all k n occur near the singularity of equation (22) (see appendix C), algebriac expressions 
can be derived (see appendix E) for the values of t n> Thus, all the quantities necessary 
for evaluating a n and a n have been determined. 

The values of the coefficients which result from incorporating the orthogonality rela- 
tion are the values of b n> Consequently, the error can readily be evaluated by comparing 

b with the excessive and deficient values of a . 
n n 


DISCUSSION OF RESULTS 


The evaluation of error due to the orthogonality assumption for the nth mode of the 
electric field requires a numerical calculation of a n , a n , and b n for a specific case. 
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Since calculation of the b n for one specific physical model had already been made 
(ref. 2) (with the assumption of orthogonality of the Bessel functions), this model was 
used herein. Only conditions that fall in regions of interest were checked. Of course, 
any region where a wave number is complex cannot be checked because the final equations 
presented herein would not be applicable. The results of reference 2 are in terms of an 
efficiency which is proportional to the nth mode of the electric field. Consequently, the 
greatest percentage error calculated in the coefficients herein would also be the greatest 
error in efficiency. 

With these conditions in mind, corrections described in this report were determined 
for the selection of points indicated in figure 2 (ref. 2). The values of a n , a n , and b n 
are presented in table I. The accuracy was sufficient to allow an error calculation to 
within at least 0. 1 percent. The results presented in table I indicate that the orthogonal- 
ity assumption of references 1 and 2 was well founded. The errors in determining the 
coefficients are insignificant for the fundamental mode, being less than 0.1 percent for 
the conditions tested. 

The only condition where there is a considerable difference is in the second and third 
modes near the first maximum in the efficiency. However, the percentage of the total 
field energy involved in the second and third modes at this condition is quite small when 
compared with the field energy of the fundamental. Thus, the quantity of primary inter- 
est in the present case, the efficiency of energy addition, is negligibly affected by the 
errors resulting from the orthogonality assumption. If the field amplitude of a specific 
higher mode were required, for diagnostic purposes, for example, the orthogonality as- 
sumption would prove unsatisfactory for some conditions. 


CONCLUDING REMARKS 

Nonorthogonal eigenfunctions are not uncommon in physical problems: they arise 
especially when wave propagation is involved. When an expansion in terms of these 
eigenfunctions is desired, the most expeditious approach could well be to use the orthog- 
onality assumption and develop a Fourier series. The resulting simplification in the 
analysis may easily justify the procedure. The accuracy of the results might be eval- 
uated by comparing the Fourier -Bessel coefficients with those determined by the numeri- 
cal methods used in this report. 

Indeed there are other examples in the literature where the orthogonality approxima- 
tion was used to obtain a Fourier expansion (ref. 7). This approach may be justifiable 
when only a slight modification of the boundary conditions will give orthogonal eigenfunc- 
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tions. The eigenfunctions for the vacuum boundary case might then be termed almost- 
orthogonal functions. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 30, 1965. 
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APPENDIX A 


SYMBOLS 


A 

nm 

defined by eq. (5) 

a 

reciprocal of square of nondi- 

f 

A 

nn 

defined by eq. (6) 


mensional Alfven velocity for 
dense plasmas (ref. 1) 

a 

deficient coefficient 

$ 

number larger than every a m 

a 

excessive coefficient 


(evaluated in appendix B) 

a m 

coefficient of mth term in 

6 nm 

Kronecker delta 

Bessel expansion of J.,(i/R) 

K 

nondimensional axial wave num- 

b n 

coefficient obtained by orthog- 


ber (ref. 1) 

onality assumption 

V 

nondimensional radial wave num 


defined by eq. (4) 


ber (ref. 1) 

C nm 

(8 nm - 1)A /A_ 

v nm ' nm' nn 

SI 

nondimensional frequency 

■*on 

nth zero of zero-order Bessel 


(ref. 1) 

function of first kind J Q 

Subscripts: 


, 2 2 l/ 2 

m 

pertaining to mth natural mode 

k n 

- 0 1 


of plasma oscillation 

R 

nondimensional radius (ref. 1) 

n 

pertaining to nth natural mode 

R 1 

nondimensional plasma radius 
(ref. 1) 


of plasma oscillation 
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APPENDIX B 


EVALUATION OF & 

In order to pick a number, which is greater than every a m for all m > N + 1, 
a set of numbers must be examined. According to reference 5 (p. 26), if 


£ fnmh 1 

m=l 

there exists a set of numbers given by 


S n > 3 

1 - £ l C „ml 

m=l 

such that the largest will satisfy the condition desired of that is, a m for 
m > N + 1. To guarantee that $ was indeed larger than any the relation 

S- 1- 1 (maxj^J) 

was used to allow for any computational errors. 

Because consists of an infinite set of elements, it is not particularly straight- 
forward to find their maximum. However, if |^J can be shown to be a monotonically de- 
creasing function of n, its maximum is just the first element of the series. Difficulty 

OO 

arises, though, in evaluating ^ l^nml because the elements of the sum contain compli- 

m=l 

cated collections of Bessel functions. As an alternative, the Bessel functions are approxi- 
mated by expanding them about the zeros of J Q . In appendix C, it is shown that the argu- 
ments will be near the zeros of J Q and that the arguments k n R^ will be near 

2 2 

(k q - O ) Rj, where k is given by equation (Cl). In addition, for higher order solu- 
tions the arguments will approach these limits more closely. This is demonstrated in the 
expansion of the Bessel function, carried out in appendixes D and E. 

In order to verify that \e^ n \ is a monotonically decreasing function, approximate ex- 
pressions for C n , A nm , and A nn must be derived. It follows from equation (D7) that 
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(Bl) 


for t n « 1 


J oK R l> 

V R ° n 1 ss -j t 

” hW 


and from equation (E4) that 



where < j j. Because 


^n R l “ 3on * ( n - 1/4 ^ 
(appendix E) equation (Bl) may be written as 


u n R l 


J o^n R l) 


« - 1 


• — V d i 

4n/ 1 


J l(^ n R l) 

Because the zeros of J Q occur near the maximum of J^, 


V ™n R l 


Equations (B2), (B3), and (B4) are the approximations that reduce C n , A nm , and A 
to algebraic expressions. Substituting these approximations and the condition n » 
into C n (eq. (4)) results in 


C n» 


2 



where F is a parameter independent of n and given by 


yR J n (^Ri) 

F = — — — — 

7 tD^ Jj(i^Rj) 


(B2) 


(B3) 


(B4) 


nn 

1 


(B5) 
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Equation (B5) reduces to 


(1+F) 


n 


. \5/2 5/2 

-;) " 


(B6) 


Similarly, A nn for n » 1 becomes 


R 1 1 

A « — - oc ± 

nn 0 

2 n 
7 t n 


Consequently, 


k n 1 

b = oc — 

n 


A 3/2 
nn n ' 


The term b n consequently has for n » 1 an absolute value that is monotonically de- 
creasing. Finally, A^ m for large argument may be written as 


A cc 

nm 


(m + „) /m - i) 1/2 m„ 3 /2 


The expression for C nm is, consequently, 


C 


nm 


(m + „) (m - l'| 1/2 m„ 1 /2 


which is a monotonically decreasing function of n. Because each term of the series 


lie 


nm I 


m=l 


is a decreasing function of n, the sum over m also is a monotonically decreasing func- 
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tion of n. Therefore, monotonically decreases for large n. To determine the var- 

iation of for small n, numerical calculations were made of jjf n \ for all n less than 
100, where the sum in included 300 terms. In the range 1 < n < 100, was also 

continuously decreasing. Consequently, the maximum value of this set occurs at a con- 
dition of n = 1. 

There remains the verification that 


I l c .J < 1 

m=l 


where the inequality must hold for all n. Expressions for A nn and A nm must be writ- 
ten that are valid in the region of small n. These are 


R?J?(j ) / 4t 2 

A „ 1 l VJ oir L n 

2 

\ ■'on 


R 


v 2 (n - -I 


4t" 


1 - 


/ l\ 2 2 

n - - n 


(B7) 


and 


r; 


A. 


nm 


2» 2 „ m ^n-iy /2 (m-iy /2 (m + n-l’ 


(B8) 


From equations (B7) and (B8) the sum may be written as 
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Because n > 1, equation (B9) reduces tc 



In order to estimate the magnitude of ^ l^nml’ *n have to be evaluated. The 

m=l 

most adverse value of t n occurs when n = 1. Because t^ never exceeds 0. 2 for the 

■1 


cases calculated, the maximum value of 11 


4*tf 

sV' 


will always be less than 1. 03, and 


consequently, 


oo 


Z 



However, from reference 8 



— = 1. 644 
2 


m 


m=l 


Therefore, 


OO 



which is less than 1. 
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APPENDIX C 


CHARACTERISTICS OF NATURAL MODES 


The natural modes are determined by the simultaneous solution of the dispersion re- 
lation (ref. 6) 


2 = « 4 0! 2 - C2 2 [2a + (k 2 - & 2 J(k 2 - fl 2 ) + (k 2 - SI 2 )' 
0 2 [<* + (/c 2 - O 2 )] - ( k 2 - S 2 2 ) 


( 22 ) 


and the boundary condition 


J 0 (^ R i) 2 2 1/2 K o[<« 2 ' ° 2 ) 1/2r J 

0 1 = ~( k z - n 2 ) Rj 0L u 


uR- 


Jl^Rj) 


K 1 [(/c 2 - J2 2 ) 1/2 rJ 


(23) 


These two equations are presented in figure 3 (fig. 2 of ref. 2) as a log log plot. 

2 

Because v obtained from equation (22) is a monotonically decreasing function with 
o 

respect to k and because the regions of possible solutions of v n (fig. 1) are located at 
larger values of v as n increases, /c n will decrease and approach the finite limit of 


K 2 = _A_ + n 2 

° l-f2 2 


(Cl) 


The right side of equation (23) consequently approaches a constant as v increases. Be- 
cause u n increases without limit and because the amplitude of is continuously de- 
creasing, ^ n Rj must approach the zeros of J Q in order for equation (23) to be satisfied. 
These characteristics of the natural modes (v n and « n ) suggest an expansion of J Q , J^, 
and u n about the zeros of J Q and an expansion of k about the value of k q given by 
equation (Cl). 
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APPENDIX D 


EXPANSION OF J Q AND Jj ABOUT ZEROS OF J Q 

It is convenient (see appendix C) to expand J 0 (^ n Rj) and J^(^ n R^) about j , which 
is the nth zero of the zero-order unmodified Bessel function of the first kind. Hence, 
because 

•ToW = J o<>on + ‘n> 

J l("n R l) = J l«on + 

and | t n J < |j Qn |, the following addition formula (ref. 9, p. 143) may be utilized: 


J o (j on + t n> = E J -m^on )J m(y 


m--°o 


W.U * 2 £ (-D-'Jm'WW 

m=l 


(Dl) 


Now, because J 0 (j on ) = 0, write from the recurrence relations (ref. 9, p. 45) that 


J 2«on> = J- J 1<W 
•'em 


k -*on 


4 w 


J 5 ( W " \~4 


r 2 7 • 3 2 3 • 5 

.2 

^on **on 


+ 1 Ww 


(D2) 
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In addition, for t small 


J (t ) = 
m tv 




2 m m! 


(D3) 


Substituting (D2) and (D3) into (Dl) and dropping all the terms of the third power or higher 
of t result in 


J o<% R l> * 



L^on^ 


(D4) 


Because the series appearing in equation (Dl) is convergent, the series from which equa- 
tion (D4) was obtained is also convergent. The addition formula for Jj(j on + t n ) may be 
written as 


oo 

^l^on + *n^ " S ^1-m^on^m^n^ 
m=-°° 


Rearranging the terms results in 


J l«on + t r)= 1, (- 1 ) m_1 [ J m-l( t n) " WV] 


^m^or) 


m=l 


(D5) 


Substituting equations (D2) and (D3) into equation (D5) results in 


J iW * 


t t 2 t 2N 

1 “_+_2 “jj (j ) 

j .2 n 2 1 on 

] ° n ^n ' 


(D6) 


Because the series appearing in equation (D5) is convergent, the series from which equa- 
tion (D6) was obtained is also convergent. Consequently, the expression l/ n ^i[j 0 ( i; n ^i)/ 
Jl^nRiilthat appears in the expression for A nrn (eq. (5)) may be written as 
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J o^n R l> _ ~ 43 o nt nyon + 2 / 
J l( I 'n R l> 4j* - 4j t + 4 1 2 „ - j* t? 


n •’on n 
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APPENDIX E 


EVALUATION OF t n 

The term t n may be found by first substituting the expansion formula of appendix D 
into the two expressions used to determine v and « n (eqs. (22) and (23)). In addition, 
similar expansion formulas for the x n variable are required before equation (23) reduces 
to a purely algebraic expression. Once this is done, the algebriac equations can be solved 
for t n . 

A characteristic of the solutions of k that is useful to evaluate higher modes is the 
property that x n is a very slowly decreasing function which has a finite lower limit. Be- 
cause x n lies close to a fixed value, it is possible to expand the right side of equation 
(23) by Taylor’s expansion. If the expansion is carried out around 


2 = _n^ 2 

° 2 

i - sr 


(El) 


(see appendix C), the result is 


‘^on^on + J 


4 ion- 4 i on t n + 4 t n-ion t n 


= -<B + 


2B - (x 2 - S2 2 )R? + B 2 ! K °- n 

° 1 J 9. ■ 


2 0 2 


where 


(3 + 2B)(x 2 - ft 2 )R 2 + B(2 + 5B + 2B 2 ) 


X 2 A 2 
o n 


/ 2 0 2 v 2 

(k o - O ) 


(E2) 


2 o ^ 2 1 

9 9 1/2 K ( k * - ST) R. 

B = (x 2 - S2 2 ) R 1 ° L 1J 


2 2 ^ 2 

Kj[(x^ - an R x 


and A n = k r - x Q . Substituting = j Qn + t n and x n = x Q + A n into equation (22) gives 
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(E3) 


gon + 

R? 


A 2 


(1 - S2 2 ) (2/c q A n + A 2 ) 


< 2 *o A n * 


The nth positive zero of J 0 (i 0n ) can readily be determined from the relation (ref. 9, 
p. 505) 



6277237 

3440640 



Equations (E2) and (E3) are two nonlinear algebraic equations in the two unknowns A n 
and t and consequently can be solved. 

A somewhat simpler expression for t will suffice if the first few natural modes 
are known. For these higher modes, an empirical relation of the form 


P 



m=l 


(E4) 


can be written. The first p values of t are obtained from 

^n R l = V + *n 

where the i^’s were obtained from reference 2. Using these values of t n in equation 
(E4) yields p simultaneous equations for the D 's. The resulting equation for t fl can 
then be used to calculate t for n’s > p. These first few modes were obtained in refer- 
ence 2 by solving equations (22) and (23) directly by a modified Newton -Raphson technique. 
However, this method (ref. 2) breaks down when used for the higher modes. The results 
from the empirical relation (E4) and from equations (E2) and (E3) indicate that both ap- 
proaches were compatible with the original equations, that is, equations (22) and (23). 
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TABLE I. - DEFICIENT, EXCESSIVE, AND APPROXIMATE COEFFICIENTS 


Ion density, 

Magnetic 

Mode, 

Deficient 

Excessive 

Orthogonal 

. / 3 

ions/ cm 

field, 

n 

coefficient, 

coefficient, 

coefficient, 


G 


*n 

*n 

b n 

10 12 

4. 2015X10 3 

i 

1. 00064 

1.00071 

1. 00068 



2 

3. 44680X10' 4 

3. 81019X10" 4 

-1. 73485X10" 3 



3 

-4. 88554 

-4. 58876 

6. 31477X10” 4 


4. 24585X10 3 

1 

5.90566X10" 1 

5. 90595X10" 1 

5.90722X10 -1 



2 

-1.75595 

-1.75580 

-1.76466 



3 

9.00958X10" 2 

9. 01084X10 -2 

9. 05804X10 -2 


4. 26X10 3 

1 

4. 91661X10' 1 

4. 91683X10" 1 

4. 91782X10" 1 



2 

-1. 51788 

-1. 51771 

-1. 52439 



3 

7. 82748X10" 2 

7. 82843X10" 2 

7.86396X10" 2 


4. 31396X10 3 

1 

1. 75898X10' 1 

1. 75904X10' 1 

1. 75931X10" 1 



2 

-5. 71984X10' 2 

-5. 71955X10' 2 

-5.73651X10" 2 



3 

2. 97105 

2. 97130 

2. 98030 

10 14 

4. 61766X10 3 

1 

1.66201X10" 1 

1. 66203X10" 1 

1. 66208X10" 1 



2 

-5. 59803X10 -2 

-5. 59779X10" 2 

-5.60098X10" 2 



3 

2. 92708 

2.92737 

2. 92893 
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Figure 1. - Schematic of Bessel functions, J 0 and Jj, 
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Figure 1. - Schematic of Bessel functions, J 0 and Jj. 
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Figure 3. - Graphic solution for natural modes of plasma 
oscillation. 
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